HYDRODYNAMICS OF RELATIVISTIC FIREBALLS 



Tsvi Piran, 1 ' 2 Amotz Shemi 3 and Ramesh Narayan 



1 



Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, USA. 



2 



Racah Institute for Physics, The Hebrew University, Jerusalem, Israel. 



3 



Department of Physics and Astronomy, Tel Aviv University, Tel Aviv, Israel. 



Summary 



We examine analytically and numerically the evolution of a relativistic fireball. We 
show that, after an early rearrangement phase, most of the matter and energy in the 
fireball is concentrated within a narrow shell. The shell propagates at nearly the speed 
of light, with a frozen radial profile and according to a simple set of scaling laws. The 
spectrum of the escaping radiation is harder at early times and softer later on. The results 
are applicable to models of 7-ray bursts. 
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1. Introduction 



The sudden release of a large quantity of gamma ray photons into a compact region 
can lead to an opaque photon-lepton "fireball" through the production of electron-positron 
pairs. The term "fireball" refers here to an opaque radiation-plasma whose initial energy 
is significantly greater than its rest mass. The formation and evolution of fireballs is of 
interest in astrophysics (Cavallo & Rees 1978), especially for the understanding of gamma- 
ray bursts at cosmological distances (Goodman 1986, Paczynski 1986, 1990, Shemi & Piran 
1990, Narayan, Paczynski & Piran 1992, Meszaros and Rees 1992a,b), or in the halo of the 
Galaxy (Piran & Shemi 1993). 

In this paper we investigate the hydrodynamics of fireballs. We begin by summarizing 
in this section several qualitative results that are already known in this problem. Because 
of the opacity due to pairs, the radiation in a fireball cannot initially escape. Instead, the 
fireball expands and cools rapidly until the temperature drops below the pair-production 
threshold and the plasma becomes transparent. In addition to radiation and e + e~ pairs, 
astrophysical fireballs may also include some baryonic matter which may be injected with 
the original radiation or may be present in an atmosphere surrounding the initial explosion. 
The electrons associated with this matter increase the opacity, delaying the escape of 
radiation. More importantly, the baryons are accelerated with the rest of the fireball and 
convert part of the radiation energy into bulk kinetic energy. 

As the fireball evolves two important transitions take place. One transition corre- 
sponds to the change from optically thick to optically thin conditions. As long as the 
total opacity (pairs + matter) is large the plasma expands adiabatically as a perfect fluid 
(Goodman 1986). However, once r drops below 1, the photons and baryons decouple from 
each other and continue their evolution independently and without interaction. The sec- 
ond transition corresponds to the switch from radiation dominated to matter dominated 
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conditions, i.e from rj > 1 to rj < 1, where rj is the ratio of the radiation energy E to the 
rest energy M: rj = E/Mc 2 (Cavallo & Rees 1978, Shemi & Piran 1990). In the early 
radiation dominated stages when rj > 1, the fluid accelerates in the process of expansion, 
reaching relativistic velocities and large Lorentz factors. The kinetic energy too increases 
proportionately. However, later when rj < 1, the fireball becomes matter dominated and 
the kinetic energy is comparable to the total initial energy. The fluid therefore coasts with 
a constant radial speed. The overall outcome of the evolution of a fireball then depends 
critically on the value of rj when r reaches unity. If rj > 1 when r = 1 most of the energy 
comes out as high energy radiation, whereas if rj < 1 at this stage most of the energy has 
already been converted into kinetic energy of the baryons. 

The opacity itself has a contribution from electron-positron pairs as well as electrons 
associated with the baryons. Initially, when the local temperature T is large, the opacity is 
dominated by e + e~ pairs (Goodman, 1986) . But this opacity, r p , decreases exponentially 
with decreasing temperature, and falls to unity when T = T p 20 KeV. The matter 
opacity, r&, on the other hand decreases only as -R -2 , where R is the radius of the fireball. 
If at the point where t p = 1, is still > 1, then the final transition to r = 1 is delayed 
and occurs at a cooler temperature. 

The initial ratio of radiation energy to mass, rji, determines in what order the above 
transitions take place. Shemi and Piran (1990) identified four regimes: (i) rji > rj pa i r = 
[3a^EiaTp /ATrrripC 4 Ri] 1 / 2 (where Ei and Ri are the initial energy and radius): In this 
regime the effect of the baryons is negligible and the evolution is of a pure photon-lepton 
fireball. When the temperature reaches T p , the pair opacity r p drops to 1 and 1. At 
this point the fireball is radiation dominated (rj > 1) and so most of the energy escapes 
as radiation, (ii) rj pa i r > Vi > Vb = (3a"T-E , i/87rm p c 2 i? 2 ) 1 / 3 : Here, in the late stages, 
the opacity is dominated by free electrons associated with the baryons. The comoving 
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temperature therefore decreases far below T p before r reaches unity. However, the fireball 
continues to be radiation dominated as in the previous case, and most of the energy still 
escapes as radiation, (iii) 775 > r\i > 1: The fireball becomes matter dominated before it 
becomes optically thin. Therefore, most of the initial energy is converted into bulk kinetic 
energy of the baryons, with a final Lorentz factor 7/ = rji + 1. (iv) r\i < 1: This is the 
Newtonian regime. The rest energy exceeds the radiation energy and the expansion never 
becomes relativistic. 

The above summary describes the qualitative features of a roughly homogeneous ex- 
panding fireball. In this paper we investigate some aspects of the evolution of inhomo- 
geneous fireballs with a non-uniform radial profile. We show in §2 that, after an initial 
rearrangement phase, the evolution is well described by an asymptotic large 7 solution. 
The radial profile of the fireball remains frozen over most of this phase, and each streamline 
follows simple scaling laws as a function of radius. In §3 we solve numerically the adiabatic 
expansion of a spherical fireball and compare the results with the asymptotic solution. We 
show that the agreement with the theoretical solution is good. Finally in §4 we summarize 
the results and discuss their implications. 



We consider a spherical fireball with an arbitrary radial distribution of radiation and 
matter. Under optically thick conditions there is strong coupling between the photons 
and baryons, and so the radiation and matter at each radius behave like a single fluid, 
moving with the same velocity. Since the radiation pressure dominates, the pressure p and 
the energy density e are related by p = e/3, and we can rewrite the standard relativistic 
conservation equations of baryon number and energy momentum (Weinberg, 1973) as: 



2. Scaling Laws and Asymptotic Solutions 
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where 7 = u*, u = u r = a/7 2 — 1, and we use units in which c = 1 and the mass of the 
particles m = 1. The mass density, n, the total energy density, e (which includes contribu- 
tions from the radiation as well as the relativistic electron-positron pairs at temperatures 
where the latter are present), and the pressure, p, are measured in the local frame of the 
fluid, but r and t are in the observer frame. 

Change variables from r, t to r, s = t — r. Equations (l)-(3) then become 
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where the derivative d /dr now refers to constant s, i.e. is calculated along a characteristic 
moving outward at the speed of light. After a short acceleration phase we expect that the 
motion of a fluid shell will become highly relativistic (7 ^> 1). If we restrict our attention 
to the evolution of the fireball from this point on, we may treat 7 _1 as a small parameter 
and set 7 ~ u, which is good to order o(7 -2 ). Then, under a wide range of conditions, 
which we discuss below, the quantities on the right-hand sides of equations (4)-(6) are 
significantly smaller than those on the left. We therefore set the terms on the right to 
zero, and obtain the following conservation laws for each fluid shell: 

r 2 n7 = constant, r 2 e 3//4 7 = constant, r 2 + —e^j 7 2 = constant. (7) 



Two regimes of behavior are then immediately apparent. In the radiation-dominated 
phase (e >> n), we have 

7 oc r, n oc r -3 , e oc r -4 , T b s ~ constant, (8) 

where T ^ s oc 7c 1 / 4 is the temperature of the radiation as seen by an observer at infinity. 
(Strictly, the radiation temperature depends on e r , the energy density of the photon field 
alone; for T <C m e c 2 , e r = e, but for T > m e c 2 , e contains an additional contribution 
from the electron position pairs, see Shemi & Piran 1990; we neglect this complication for 
simplicity). The scalings of n and e given in (8) correspond to those of a fluid expanding 
uniformly in the comoving frame. Indeed, all four scalings in eq (8) were derived for a 
homogeneous radiation dominated fireball by Shemi & Piran (1991, see also Goodman 
1986) by noting the analogy with an expanding universe. What we have shown here is 
that the same relations are valid for each individual radial shell in the fireball even in the 
more general inhomogeneous case. In fact, these scaling laws also apply to Paczyhski's 
(1986) solution for a steady state relativistic wind. When we neglect the right hand sides 
of eqs. (4)-(6) the problem becomes effectively only r dependent. 

Although the fluid is approximately homogeneous in its own frame, because of Lorentz 
contraction it appears as a narrow shell in the observer frame, with a radial width given 
by Ar ~ r/7 ~ constant ~ Ri, where Ri is the initial radius of the fireball. We can now 
go back to eqs (4)-(6) and set d/ds ~ 7/r. We then find that the terms we neglected on 
the right hand sides of these equations are smaller than the terms on the left by a factor 
~ I/7. Therefore, the conservation laws (7) and the scalings (8) are valid so long as the 
radiation-dominated fireball expands ultra-relativistically with large 7. The only possible 
exception is in the very outermost layers of the fireball where the pressure gradient may 
be extremely steep and d/ds may be ^> 7/r. Ignoring this minor deviation, we interpret 
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eq (7) and the constancy of the radial width Ar in the observer frame to mean that the 
fireball behaves like a pulse of energy with a frozen radial profile, accelerating outward at 
almost the speed of light. 

In the alternate matter-dominated regime (e < n), we obtain from eq (7) the following 
different set of scalings, 

7 — > constant, nocr -2 , e oc r~ 8 / 3 , T Q b s oc r~ 2//3 . (9) 

The modified scalings of n and e arise because the fireball now moves with a constant 
radial width in the comoving frame. (The steeper fall-off of e with r is because of the 
work done by the radiation through tangential expansion.) Moreover, since e <C n, the 
radiation has no important dynamical effect on the motion and produces no significant 
radial acceleration. Therefore, 7 remains constant on streamlines and the fluid coasts with 
a constant asymptotic radial velocity. Of course, since each shell moves with a velocity 
that is slightly less than c and that is different from one shell to the next, the frozen 
pulse approximation on which eq (7) is based must ultimately break down at some large 
radius. We consider this question below, but first continue with our investigation of the 
approximate relations in eq (7). 

A scaling solution that is valid in both the radiation-dominated and matter-dominated 
regimes, as well as in the transition zone in between, can be obtained by combining the 
conserved quantities in eq (7) appropriately. Let to be the time and ro be the radius 
at which a fluid shell in the fireball first becomes ultra-relativistic, with 7 > few. Label 
various properties of the shell at this time by a subscript 0, e.g. 70, no, eo, and define 
77 = e/n, t]o = eo/no. Defining the auxiliary quantity D, where 

I^7o + j37o__ A, (10) 
D 7 4?7o7 4 ?7o 
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These are parametric relations which give r, n, e, and r\ of each fluid shell at any time 
in terms of the 7 of the shell at that time. The relation for r in terms of 7 is a cubic 
equation. This can in principle be inverted to yield 7(r), and thereby n, e, and rj may also 
be expressed in terms of r. 

The parametric solution (11) describes both the radiation-dominated and matter- 
dominated phases of the fireball within the frozen pulse approximation. For 7 <C ?7o7o, 
the first term in eq (10) dominates and we find D oc r, 7 oc r, recovering the radiation- 
dominated scalings of eq (8). This regime extends out to a radius r ~ ?7o r o- At larger 
radii, the first and last terms in (10) become comparable and 7 tends to its asymptotic 
value of 7/ = (4?7o/3 + 1)70- This is the matter dominated regime. (The transition occurs 
when 4e/3 = n, which happens when 7 = 7//2.) In this regime, D oc r 2 / 3 , leading to the 
scalings in eq (9). 

Ultimately, all the energy in the fireball is concentrated in the kinetic energy of the 
matter, and this determines the value of 7/. Interestingly, if we write 7/ in terms of the 
initial parameters of the fireball at time t = 0, we find 7/ = rji + 1, whereas, when we 
write it in terms of i] , 70, we have the additional factor of 4/3 as written above. Both 
formulae represent energy conservation, but the component T u of the energy momentum 
tensor behaves differently in the two cases. At time t = 0, the fluid is at rest and the 
radiation energy density is merely e, whereas at t = t , the fluid is already moving highly 
relativistically and there is an additional contribution to the energy from the moving 
pressure, T u = + (l/3)ujje ~ (4/3)7 2 e. 

Let us now return to a consideration of very late times in the matter-dominated phase 
when the frozen pulse approximation begins to break down. We have already seen that 



in this phase the radiation density e is much smaller than the matter density n, and also 
that 7 tends to a constant value 7/ for each shell. We may therefore neglect the term 
— (l/3)(<9e/<9r) in eq (3) and treat 7 and u in eqs (l)-(3) as constants. We then find that 
the flow moves strictly along the characteristic, 13 ft — r = constant, so that each fluid 
shell coasts at a constant radial speed, /?/ = M//7/. Let us label the baryonic shells in 
the fireball by a Lagrangian coordinate R, moving with a fixed Lorentz factor 'jf(R), and 
let t c and r c represent the time and radius at which the coasting phase begins, which 
corresponds essentially to the point at which the fluid makes the transition from being 
radiation dominated to matter dominated. We then find 
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The separation between two neighboring shells separated by a Lagrangian distance AR 
varies during the coasting phase as 
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Thus the width of the pulse at time t is Ar(t) pa Ar c + A^f(t — t c )/7/ ~ Ri + (t — *c)/7/> 
where Ar c ~ Ri is the width of the fireball when it begins coasting, 7/ is the mean 7/ in 
the pulse, and A7/ ~ 7/ is the spread of 7/ across the pulse. From this result we see that 
there are two separate regimes in the fireball evolution even within the matter dominated 
coasting phase. So long as t — t c < ^jRi, we have a frozen-coasting phase in which Ar 
is approximately constant and the frozen pulse approximation is valid. In this regime the 
scalings in eq (9) are satisfied. However, when t — t c > 'jfRi, the fireball switches to an 
expanding-coasting phase where Ar oc t — t c and the pulse width grows linearly with time. 
In this regime the scaling of n reverts to n oc r -3 , and, if the radiation is still coupled to 
the matter, e oc r -4 . 



Independently of the above considerations, at some point during the expansion, the 
fireball will become optically thin and the radiation will decouple. From this stage on the 
radiation and the baryons no longer move with the same velocity and the radiation pressure 
vanishes, leading to a breakdown of equations (2) and (3). The radiation will now coast 
with a speed exactly equal to c and with a constant radial width. The radiation energy 
density will clearly scale as e oc r~ 2 . The baryon shells on the other hand will coast with 
their own individual velocities. If the fireball is already in the matter dominated coasting 
phase there will be no change in the propagation of the baryons. However, if the fireball 
is in the radiation dominated phase when it becomes optically thin, then the baryons will 
switch immediately to a coasting phase. 

3. Numerical Results 

We have developed a spherically symmetric relativistic code that follows the evolution 
of a mass-loaded fireball from an initial configuration at rest via the acceleration phase 
into the asymptotic frozen pulse regime. The code is Eulerian and employs a second 
order conserved scheme (Bowers and Wilson, 1991) modified to take care of the extreme 
relativistic Lorentz factors encountered in this problem. Our scheme is quite different 
from the one used by Vitello and Salvati (1976), who studied a similar problem. The code 
has passed several standard tests, including the Richardson test, i.e. the results converge 
satisfactorily as the grid size is decreased. 

The initial profiles for the cases that we present here are: 

e(r,f = 0) = J^s-y ™M) = ^, 7M) = 1, (14) 

where we choose Ri, the initial width of the fireball, to be unity. The initial radiation 
density eo is in arbitrary units, and we assume a constant ratio r\i of energy to mass. The 
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energy density falls off sufficiently rapidly with radius in the initial profile that the external 
density is negligible compared to the interior density. However, we cannot set the exterior 
density exactly to zero since this leads to numerical problems. We have explored different 
initial conditions, and find different pulse shapes, but the overall qualitative behavior is 
generally similar to the results described below. 

Fig. l.a shows a sequence of profiles in the observer frame of the energy density, 7e, 
and the mass density, 777-, for a simulation with r\i = 50. Three phases of evolution are 
apparent, viz. an initial acceleration/rearrangement phase, a short radiation dominated 
phase, and a final matter dominated phase. Conservation of energy requires that, asymp- 
totically, the average Lorentz factor of the expanding fluid should be 7/ = r\i + 1. We 
compute 7 at each stage by means of 



where the second expression is preferable for numerical accuracy. The average Lorentz 
factor in this simulation does approach the expected asymptotic value, but it does not 
quite reach it because we did not let the computation continue for a long enough time. 

Early in the evolution, the pulse rearranges itself during a brief acceleration phase. 
Although with our choice of initial data the ratio between the energy and mass density, 
?7, begins with a constant value throughout the fireball (eq 14), r\ changes during this 
phase and no longer remains constant within the pulse. Generally, r\ ends up being smaller 
in the inner parts of the fireball and larger on the outside. This can be seen in Fig. 2 
which shows n, e and 7 at a fixed moment of time after the end of the early phase of 
rearrangement. After this phase, the shape of the pulse is frozen and the fireball evolves 
through a radiation-dominated phase to a matter-dominated state. The transition from 
radiation dominated to matter dominated can be clearly seen in Fig. l.a where the matter 
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jT u r 2 dr 



jT u r 2 dr 
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density is initially lower than the energy density but becomes larger by the end of the 
computation. 

The profile of the Lorentz factor 7 at various times in the pulse is shown in Fig. l.b. 
We see that 7 varies significantly across the pulse. Whereas the mean 7/ cannot exceed 
the asymptotic value of rji + 1, the maximum Lorentz factor within the pulse is larger and 
in fact increases throughout the evolution. This happens as the outermost layers of the 
fireball keep accelerating, as can be seen in Fig. 2. The Lorentz factor therefore peaks 
ahead of the energy density in a low density region, and a small fraction of the material is 
accelerated to these high 7 values. The peak in 7 leads to the highest observed temperature 
coming from the front of the fireball and to lower temperatures coming from the interior 
(Fig. l.b). 

In Fig. 2 we compare the calculated pulse at t = 153.7 to a pulse extrapolated using 
eqs. (10)-(11) from t = 14. The agreement is very good considering that the energy and 
matter density have decreased by four and three orders of magnitude, respectively, while 
the maximal Lorentz factor has increased by a factor of six between these two moments. 
The agreement is excellent in the inner region of the pulse but is less satisfactory on 
the outside. This is partly due to continued acceleration of the outer regions where the 
pressure gradient is steepest, and partly because of decreased numerical accuracy in the 
regions where the density is very low. 

Fig. l.b shows profiles of T Q b s , the temperature that would be observed at infinity if 
the radiation could escape. At early times, there is a drop in T b s due to the broadening 
of the pulse as a result of internal rearrangement. Then, during the subsequent radiation 
dominated phase, T oc 7 _1 (eq 8), and the observed temperature T obs of each shell remains 
a constant. The overall spectrum of a fireball that becomes optically thin during this 
phase is a blending of thermal spectra with different temperatures and different blueshifts, 
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and will be slightly broader than a single temperature thermal spectrum (see Goodman, 
1986). The spectrum does not change during the radiation dominated phase, apart from 
a minor effect due to the addition of a harder component from the acceleration of the 
outermost layers of the fireball. When a given shell enters the matter dominated phase 
T b s begins to decrease, since T continues to decrease but now without a compensating 
increase in 7 (eq 9). The result is a softer spectrum than that observed during the radiation 
dominated phase. The emitted spectrum depends now on the moment that each shell 
becomes optically thin. Since different shells become optically thin at different values of 7 
we expect the spectrum to be borader than that emitted during the radiation dominated 
phase. 

The evolution of a relativistic radiation fireball that we have described here and in 
§2 is remarkably different from that of a Newtonian fireball with ^ C 1. Fig. 3 shows 
the energy density and the mass density for a pulse with r\i = 0.001. None of the features 
described earlier appear. Relativistic velocities are never reached and the shape of the 
pulse is not frozen. Instead we observe an expanding, almost homogeneous sphere, rather 
than an expanding shell of matter and radiation, and the expansion velocity of most of the 
fireball is roughly the Newtonian velocity y/2rjl = 0.045. A negligible fraction of the matter 
on the surface is accelerated to higher speeds. Interestingly, the Newtonian fireball bears 
a strong qualitative resemblance to the relativistic fireball in the local frame. Therefore, 
the differences between the two cases arise mainly because of the transformation to the 
observer frame. In the Newtonian case there is no difference between the observer frame 
and the matter local frame, but in the relativistic case Lorentz contraction leads to a 
drastic change in the appearance of the fireball. Therefore, in the former case the fireball 
appears to fill the entire sphere of radius r whereas in the latter case the observer sees a 
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narrow pulse whose width remains of the same order as the original width, leading to a 
time scale ~ Ri/c. 

4. Conclusions 

We have shown in this paper that fireballs with a large initial ratio rji of radiation 
energy to rest mass energy show certain common global features during their expansion and 
evolution. After a short initial acceleration phase, the fluid reaches relativistic velocities, 
and the energy and mass become concentrated in a radial pulse whose shape remains frozen 
in the subsequent expansion. The motion is then described by an asymptotic solution 
(eqs 10, 11, §2), which gives for each individual shell scaling laws similar to those of a 
homogeneous sphere. 

The expanding fireball has two basic phases: a radiation dominated phase and a matter 
dominated phase. Initially, during the radiation dominated phase the fluid accelerates with 
7 oc r for each Lagrangian shell. The fireball is roughly homogeneous in its local rest frame 
but due to the Lorentz contraction its width in the observer frame is Ar ~ Ri, the initial 
size of the fireball. When the mean Lorentz factor of the fireball becomes 7 ~ (rji + l)/2 a 
transition takes place to the matter dominated phase. Ultimately, all the energy becomes 
concentrated in the kinetic energy of the matter, and the matter coasts asymptotically with 
a final Lorentz factor 7/ !=s (rji + 1). The matter dominated phase is itself further divided 
into two sub-phases. At first, there is a frozen-coasting phase in which the fireball expands 
as a shell of fixed radial width in its own local frame, with a width ~ 7/-Ri ~ rjiRi. Because 
of Lorentz contraction the pulse appears to an observer with a width Ar Ri. Eventually, 
the spread in 7/ as a function of radius within the fireball results in a spreading of the pulse 
and the fireball enters the coasting-expanding phase. In this final phase, Ar ~ Rit/jp 
and the observed pulse width increases linearly with time. 
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The fireball can become optically thin in any of the above phases. Once this happen 
the system ceases to behave like a fluid, and the radiation moves as a pulse with a constant 
width, while the baryons enter a coasting phase like the one described above. 

We have verified many of these theoretical results by means of numerical simulations of 
spherically symmetric relativistic fireballs (§3). In particular, we confirm that the asymp- 
totic solution with frozen pulse shape is reproduced to a good approximation. This is a 
very useful result since it implies that in future it is not necessary to carry out numerical 
simulations to very late times. As soon as the Lorentz factor of the expanding fluid reaches 
a moderately large value, say 7 ~ 10, we can use the theoretical results to extrapolate the 
pulse. This will provide a huge saving in computation, particularly in cases where rji ^> 1 
and the asymptotic 7 is very large. 

An important aspect of fireball evolution that can be studied only by numerical sim- 
ulations is the early stages of rearrangement. During this phase the fireball is still only 
mildly relativistic and it internally modifies the profiles of energy and matter density. From 
a number of simulations with different initial conditions, we find that the ratio of energy 
to matter density, 77, usually ends up with a lower value in the interior of the fireball and 
with large values on the outside. The Lorentz factor 7 also invariably increases from the 
inside out. These modified profiles enter the frozen pulse phase and then do not change any 
further. Consequently, it appears to be a generic feature that any radiation that escapes 
from the fireball will be hot on the outside and cooler on the inside. In other words, the 
observed radiation pulse will tend always to have a spectral profile showing a characteristic 
hard-to-soft transition as a function of time. This effect will be enhanced if the early radi- 
ation from the outside is emitted in the radiation-dominated phase and the later radiation 
from the interior is released from matter-dominated layers. The hard-to-soft signature 
will be even stronger in this case. Even if the radiation is not obtained directly from the 
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fireball, but through shock re-radiation as the fireball interacts with external matter, this 
feature should still be present. 

The discussion in this paper has been restricted to general issues related to the evo- 
lution of relativistic fireballs. The most immediate application of these results is to cos- 
mological and Galactic halo models of gamma-ray bursts. Although these models differ in 
their explanations of the origin of the gamma-rays, all of them involve a stage where the 
initially injected energy goes through a fireball phase. Therefore, the scaling laws that we 
have written down for the matter and energy density, the temperature, and the Lorentz 
factor 7 will be relevant. Also, the hard-to-soft spectral evolution described above should 
be observed in each sub-burst and possibly across the whole burst as well. In fact, this pre- 
diction is probably valid regardless of how the final observed radiation is produced, whether 
it be through direct emission from the fireball when it becomes optically thin (Goodman 
1986, Paczyhski 1986, Shemi & Piran 1990) or through shock re-emission (Meszaros & 
Rees 1992a,b). 

This work was supported by NASA grant NAG 5-1904 to Harvard University and by 
a BSF grant to the Hebrew University.. 
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Figure Captions 

Fig. 1: (a) Energy density, e"f (solid lines), and mass density, 717 (dotted lines), in the observer 
frame for a numerical simulation where the initial energy to mass ratio rji = 50. b. 
Lorentz factor, 7 (solid lines), and observed temperature, T b s = 7T = 7c 1 / 4 (dotted 
lines). The temperature scale is arbitrary. 
Fig 2: Calculated (solid lines) and extrapolated (dotted lines) n, e and 7 profiles at t = 153.7 
(the end of the computation). The extrapolation is from t = 14 (shortly after the end 
of the rearrangement phase) using eqs. (10)-(11). The agreement is best in the trailing 
edge of the pulse in the interior of the fireball and is less satisfactory in the leading 
edge. This is because of the combined effects of the steep pressure gradient and the 
loss of numerical accuracy at lower densities. 

Fig. 3: Energy density, e (solid lines), and mass density, n (dotted lines), for a Newtonian or 
non-relativistic fireball with rji = 0.001. Note that, while the final time is the same 
as in Fig. 1, the pulse has propagated a much shorter distance. There is no hint of a 
shell structure or a frozen pulse shape here, in contrast to the relativistic case shown 
in Fig. 1. 
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